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Context. Understanding dark energy and measuring the topology of the Universe are two of the biggest open questions 
in physical cosmology. It was previously shown that multiple connectedness, via the twin paradox of special relativity, 
provides a novel physical justification for an assumption of the standard FLRW model: it implies a favoured space-time 
splitting (comoving coordinates). 

Aims. Could cosmic topology also imply dark energy? 

Methods. We use a weak field (Newtonian) approximation of gravity and consider the gravitational effect from distant, 
multiple copies of a large, collapsed (virialised) object today (i.e. a massive galaxy cluster), taking into account the 
finite propagation speed of gravity, in a flat, multiply connected universe, and assume that due to a prior epoch of fast 
expansion (e.g. inflation), the gravitational effect of the distant copies is felt locally, from beyond the naively calculated 
horizon. 



■ Results. We find that for a universe with a x spatial section, the residual Newtonian gravitational force (to 



o 



first order) provides an anisotropic effect that repels test particles from the cluster in the compact direction, in a way 
algebraically similar to that of dark energy. For a typical test object at comoving distance x from the nearest dense 
nodes of the cosmic web of density perturbations, the pressure-to-density ratio w of the equation of state in an FLRW 
^ universe, is ui ~ ~(x/^)''i where L is the size of the fundamental domain, i.e. of the Universe. Clearly, <g 1. For a 

I ' spatial section of exactly equal fundamental lengths, the effect cancels to zero. For a spatial section of unequal 

O ' fundamental lengths, the acceleration effect is anisotropic in the sense that it will tend to equalise the three fundamental 

^ , lengths. 

C/3 ■ Conclusions. Provided that at least a modest amount of inflation occurred in the early Universe, and given some other 

conditions, multiple connectedness does generate an effect similar to that of dark energy, but the amplitude of the effect 
at the present epoch is too small to explain the observed dark energy density and its anisotropy makes it an unrealistic 
candidate for the observed dark energy. 
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1. Introduction S-manifold of comoving space, rather than by an "infinite" 

flat space (e.g. Lmiiinet et al. 2003; Roukema et al. 2004; 

Two of the biggest open questions in physical cosmology are Zurich et al. 2005a, 2005b; Gundermann 2005). 

(1) the interpretation of the cosmological constant (or some mi i 

other form of dark energy), and (2) the global shape of the , ^hese anabases may or may not be supported by future 
Universe, including both curvature and topology. The for- observations. One of the most independent tests of the PDS 
mer is empirically known to exist from many different, more hypotheses will be high precision estimates of the total dcn- 
or less independent, observational approaches (the cosmic ^^^^ parameter ritot- The hypotheses of these authors geo- 
concordancc model, with matter density « 0.3, cosmo- metrically require Utot to be strictly greater than umty and 
logical constant/dark energy « 0.7, Hubble constant ptot should he m the predicted range 1.01 < l^tot ^ 1-02 
Ho - 70 km/s/Mpc, baryon density « 0.05, (e.g. Lahav (different authors using different methods make slightly dif- 
& Liddle 2004) ), and observations have also suggested an ^^-ef predictions)^ Future observations which yield, for ex- 
answer to the latter, but the evidence is far from being ^"^P^^' "tot = 1.001 ± 0.001, would rule out these hypothe- 
conclusive, and, in principle, a final answer may be totally ^^^^ ^ a significance level of 9a. 

beyond the reach of any observations. Independently of the present hypotheses, the question 

The analysis of empirical data that suggests the shape of the physical consequences of cosmic topology will remain, 

of the Universe consists of recent analyses of the cosmic mi- We know that locally, geometry and density are directly 

crowave background observations by the WMAP satellite, related, via general relativity. But we have very few hints 

These analyses found the temperature fluctuation map to as to how global geometry should relate to other physical 

be better modelled by a multiply-connected model of the properties. General relativity is a local theory, not a global 

Universe, for a Poincare dodecahedral space (PDS) as the one. 
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One known effect is that cosmic topology defines the 
comoving reference frame. 

It was shown by Uzan et al. (2002) and Barrow & Levin 
(2001) that resolving the twin paradox of special relativity 
in a multiply-connected universe is different than in simply 
connected Minkowski space. Moreover, multiple connected- 
ness implies a favoured space-time splitting. The authors 
show that this must correspond to the comoving reference 
frame. 

In other words, this provides a global geometrical mo- 
tivation for Weyl's postulate, i.e. the postulate that "the 
world lines of the galaxies form a three-bundle of non- 
intersecting geodesies orthonormal to a series of space-like 
hypersurfaces" . 

Could it be possible that cosmic topology could also 
provide a simple explanation for the observed density of 
dark energy? 

If there is such a mechanism, then this would most nat- 
urally be able to explain the fact why dark energy starts 
to dominate at the present epoch if this mechanism were 
linked to a property of the present cosmological epoch. 

One property of the present epoch is that density per- 
turbations have collapsed due to gravity and formed viri- 
alised structures, of which the most dense at the present 
epoch are galaxy clusters, the largest which typically are of 
mass M ~ 1O^^M0, which are formed from most of the ma- 
terial within a region of linear comoving size 10 — 20 Mpc. 

What property of virialised structures is different be- 
tween the simply connected case and the multiply con- 
nected case? 

The difference is that in the covering space (apparent 
space), the spatial distribution of distant objects is homoge- 
neous and uncorrelated with the local distribution if space 
is simply connected, while if space is multiply connected, 
then the spatial distribution of distant objects in the cov- 
ering space is not random. 

Given the non-random distribution of multiple topologi- 
cal images of objects such as clusters in the covering space, 
and the possibility that there was causal contact between 
them at some time in the past (e.g. due to some moderate 
amount of inflation), is it possible that the residual grav- 
itational effect due to distant, multiple topological images 
of a single object yields a gravitational effect different from 
that normally expected from the assumption of large-scale 
unperturbed, homogeneity? 

In Sect. [21 the assumptions made in considering this 
gravitational effect, in the covering space, between multi- 
ple copies of a "small" , massive object, perturbed from its 
initial position, are listed. 

In Sect. [3l the effect is calculated. 

In Sect, m the results of this calculation are discussed, 
and in Sect. [5] we conclude. 

For a short, concise review of the terminology, geometry 
and relativistic context of cosmic topology, see Roukema 
(2000) (this is now slightly outdated, but is sufficient for 
beginners). For in-depth review papers see, e.g. Lachicze- 
Rey & Luminet (1995); Luminet (1998); Starkman (1998); 
Luminet & Roukema (1999); workshop proceedings are in 
Starkman (1998) and following articles, and Blanloeil & 
Roukema (2000). For comparison and classification of dif- 
ferent observational strategies, see e.g. Uzan et al. (1999); 
Luminet & Roukema (1999); Roukema (2002); Rebougas 
& Gomero (2004). What might be considered as one of 
the most striking theoretical results, in the sense of provid- 



ing a direct link between the FLRW model and multiple- 
connectedness, is the implication of a favoured reference 
frame, which must coincide with the comoving reference 
frame: see Barrow & Levin (2001) and Uzan et al. (2002). 

Note that some previous work has been done on possible 
links between topology and dark energy, via the Casimir ef- 
fect. See e.g. Ahmadi & Nouri-Zonoz (2005) and references 
therein. Lachicze-Rey (1999) estimated that for a topol- 
ogy, the effect would be at least 50 orders of magnitude too 
weak. 

2. Assumptions and calculational choices 

Initially, in Sect. 12. ll we consider the se//-gravity of a clus- 
ter to itself, provided that the cluster has been perturbed 
from its initial, comoving position. In Sect. 12.21 the more 
realistic case of essentially stationary (in comoving coordi- 
nates) clusters distributed in the nodes of the cosmic web 
and the effect that their multiple images have on small "test 
objects" nearby, is considered. 

2.1. Self-gravity of a large cluster 
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Fig. 1. Schematic diagram showing the potential to a single, 
massive object in a multiply connected universe of comoving 
size L, which is perturbed to the right (of the figure) by a small 
physical distance x. The information about the changes in the 
potential generated by this object travels outward from the ob- 
ject at a finite speed, cq-\y. If L is large (about a Hubble length 
in size), then at or close to the position of the object, the poten- 
tial due to the (distant) topological image remains that of the 
unperturbed topological image, until a Hubble time has passed. 
Anthropomorphically, we could say that during about a Hubble 
time following the perturbing event, the object "does not be- 
lieve" that its topological image has been perturbed. Solid thick 
curves show the perturbed potentials, dashed thick curves show 
the original potentials, wavy thick curves symbolically show 
the information about the changed position of the object be- 
ing transmitted by gravitational waves. 

The following assumptions and choices are made: 

(1) Newtonian approximation of gravity 

(2) a flat covering space, 

(3) calculations are made in the covering space 

(4) a region at least a few times the size of the injectivity 
diameter [assumed to be on a size scale approximately 
similar to that of the diameter of the surface of last 
scattering (SLS)] has been in causal contact (e.g. due 
to inflation) 
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(5) the gravitational potential induced from density per- 
turbations at large. supra-SLS distance scales mostly 
cancels out due to homogeneity, but since multiple topo- 
logical images are non-randomly distributed, their total 
contributions to the potential may not fully cancel and 
must be calculated explicitly 

(6) up to a few tens of Mpc around a big cluster, the only 
non-uniform, long-distance contributions to the local 
potential are from its multiple topological images 

(7) time lapse assumption: the potential from these distant 
images corresponds to the state of these images as they 
were at cosmological time ti, about a Hubble time in 
the past {ti < ^recombination ^ ^o); at that cosmolog- 
ical time (possibly pre-inflation, see (4) above), these 
images consisted of not-yet-collapsed density perturba- 
tions, which had not yet had time to be significantly dis- 
placed towards other dense regions — both this assump- 
tion and the following (8) make the physically standard 
assumption that the speed of transmission of gravita- 
tional waves CQW is equal to that of the special and 
general relativistic space-time constants and the elec- 
tromagnetic transmission speed: cqw = '^ST = cqr = 
cem; see Ellis & Uzan (2005). 

(8) time lapse corollary: if the cluster is spatially offset 
(perturbed) from its original position, then the long- 
range contributions to the gravitational potential are 
only felt locally after at least a Hubble time: on a time 
scale much less than a Hubble time, local calculations 
can validly assume that only the cluster moves, not its 
topological images, since they have not yet received the 
information that the topological images have moved (see 
Fig.[T]). 

(1), (2): Assumptions (1), (2) are identical to those 
made for most cosmological A'^-body simulations of galaxy 
formation in the FLRW model, e.g. Roukema et al. (1997), 
Bagla (2005) and references therein. In other words, the 
"volume effect" and the "backreaction effect" are consid- 
ered negligible - see Buchert & Carfora (2003) for the basic 
equations. (The volume and backreaction effects correct for 
the fact that the Universe does not have an exactly homo- 
geneous FLRW metric, it is only approximately FLRW.) 

(3) : Particle-Mesh (PM) cosmological A-body codes of 
galaxy formation (and also those which combine PM on 
a large scale with alternatives on smaller scales) almost 
always assume a spatial 3-manifold which is the three-torus, 
T'^ . so calculations are made in the fundamental domain — 
there is no need to use the covering space, so point (3) is 
not needed for these type of A^-body codes. 

Most direct A^-body codes and tree codes (TC) also as- 
sume a 3-torus model, but make calculations in the covering 
space rather than in the fundamental domain. This is the 
choice (3) made here: calculations in the covering space are 
(usually) geometrically simpler in the covering space than 
in the fundamental domain. 

(4) : If there is no causal contact beyond the SLS, then no 
effect from the topologically lensed images can occur. Here, 
the case that the causal radius is much greater than the 
radius of the SLS, e.g. due to an earlier, moderate amount 
of inflation, is considered. Linde (2004) has recently argued 
that for zero or negative curvature, multiply-connected uni- 
verses are more likely than simply-connected universes, and 
that these multiply-connected universes are expected to 
have undergone a moderate amount of inflation. 



(5), (6): While the assumption that most large-distance 
effects on the potential should approximately cancel each 
other is likely to be a good approximation, the possible non- 
cancelling of the potential due to large-distance multiple 
images of a single cluster is likely to be of a similar order of 
magnitude to the effects which we are assuming to cancel 
under assumptions (5) and (6). 

Nevertheless, we are interested in investigating whether 
any long-distance gravitational effect, in addition to contri- 
butions from local inhomogeneities, occurs due to topolog- 
ical imaging. 

If the result were a large effect, then we would have to 
verify that it is fully self-consistent with the effects from 
"randomly" distributed objects. 

To some degree, we could expect that the effect has al- 
ready been partially modelled in PM and TC cosmological 
A^-body simulations, without making the assumptions (5) 
and (6). However, these simulations generally make "re- 
alistic" assumptions for long-distance gravitational effects, 
which mean some combination of assuming long-distance 
homogeneity and assuming an infinite speed cq\y of the 
transmission of gravity, i.e. ignoring assumptions (7) and 
(8). 

(7) : The time lapse between transmission of information 
from a gravitational source (density perturbation) and its 
arrival at a "target" point in comoving space is normally ig- 
nored in A^-body simulations of galaxy formation: gravity is 
assumed to be transmitted instantaneously. This is usually 
a reasonable approximation, since the gravitational effect 
from distance sources can generally be approximated (e.g. 
as in the top-down tree code simulations) by considering 
the mass in a large, distant cube of space, which subtends 
a small angle at the "target" point, as a single, very massive 
point object. 

However, for topological gravity effects, we do take into 
account this time lapse. 

(8) : The initial density perturbations from which a clus- 
ter formed are "not aware" of the fact that the cluster later 
on collapses gravitationally and moves towards a neigh- 
bouring potential well. 

Assumption (8) gives us the first of the two effects which 
we will calculate. 

First consider a static cluster, which does not move. 
This is separated from one of its "adjacent" topological im- 
ages by a "generator" . In general, this is an isomorphism 
in the covering space. In the simplest 3-torus case, it is a 
translation, and can be thought of as a vector in Euclidean 
3-space (the covering space). It is also separated from one 
of its adjacent topological images in the opposite direction, 
i.e. by the generator —gx- 

Now consider a cluster which has moved a small distance 
from its "initial" location in comoving space a short time 
ago. The contributions to the gravitational potential, near 
to this cluster, from the other topological images — which 
are distant — remain the same as they were before the 
perturbation occurred. This is shown in Fig. [TJ 

If we consider the potentials close to the object to be 
the potentials relevant for making local calculations about 
acceleration (or worldlines), then similarly to the way lumi- 
nosity distance is defined in terms of the observational 
flux / and intrinsic luminosity L, f = L/{AiTd\), we can de- 
fine the effective comoving distance rfj^g- (for a closer anal- 
ogy with luminosity distance, we could also call this the 
"gravity distance" ) in terms of the intrinsic mass m and 
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Fig. 2. Self-gravity of a large cluster. A flat, toroidal Universe 
model of comoving side length L, filled with what is assumed 
to be a homogeneous density field except for one massive, col- 
lapsed object (e.g. cluster of galaxies) shown as a black spot 
surrounded by an empty sphere from which the matter forming 
it was originally distributed. This object is slightly perturbed, 
by physical distance x, from its original position towards one of 
its adjacent virtual copies in the comoving ("apparent") space. 
Since the potential in the "central" copy of the fundamental 
domain is determined by the two adjacent copies of the funda- 
mental domain, the object "perceives" the adjacent topological 
images in their original positions. 



the component of locally felt acceleration x due to the dis- 
tant object (the "observed" gravitational acceleration), i.e. 
dg-g- is the distance satisfying 



Gm 
cS 



(1) 



In other words, d^Q is the comoving distance to a distant 
object implied by the local shape of its potential, taking into 
account the standard value of cq^ w 3 X lO^m/s rather 
than the naive Newtonian approximation cqw = ^so- 

Using this concept of distance, a cluster which has 
moved from its "initial" location in comoving space is nearer 
to one (or two or three, depending on the direction of mo- 
tion) of its topological images, since in standard physics, 
cgw is finite. 

Hereafter, we use the effective comoving distance unless 
otherwise stated. 

Let us consider the component of its motion towards one 
of its three adjacent images, so that the cluster is displaced 
by distance x (in "physical" coordinate units) from its "ini- 
tial" location in the direction of generator of comoving 
length L, as shown in Fig. [21 Since we assume Newtonian 
gravity at the present epoch, L is also the relevant distance 
in proper units. 

Given points (1) to (8), we now have a cluster which 
feels unequal gravitational pulls from a pair of its closest 
topological images: it is slightly less than L from one image, 
and slightly more than L from the opposite image. The net 
result should be a gravitational pull towards the former. 

This self-gravity effect should be absent in any A^-body 
simulation which assumes gravity is transmitted instanta- 
neously, since the effects of the two adjacent topological 
images will always perfectly cancel if cq-w = c» is assumed. 

2.2. Effect relative to the cosmic web 

We know from linear perturbation theory of the collapse 
of linear overdensities in an FLRW universe, such as the 
Zel'dovich approximation and A'^-body simulations, that 
matter generally "falls" from low density regions (voids) 
into filaments and streams along filaments towards knots 
where the filaments join together into what correspond to- 
day to massive galaxy clusters. 



In general, the less massive objects move faster than 
the higher mass objects, due to conservation of momentum 
(Newtonian): relatively low mass objects fall (in comov- 
ing coordinates) towards the massive galaxy clusters at the 
"knots" of the cosmic web, while the most massive clusters 
have relatively little peculiar velocities with respect to the 
comoving frame. 



L-x 
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Fig. 3. Effect relative to the cosmic web. As in Fig. [2] except 
that the big black spot representing a massive, collapsed object 
is no longer perturbed, and instead, a small test object located 
physical distance x from the cluster, along the line separating 
the two clusters, is shown. 



In this case, we can consider a massive cluster which 
is approximately stationary and x to be the distance re- 
maining between a test particle (of negligible mass) and the 
cluster. The test particle is so far mostly comoving with the 
comoving reference frame, i.e. its peculiar velocity (velocity 
relative to the comoving frame), at which it falls towards 
the cluster is "small" . 

The same assumptions and choices, (1) to (7), are made 
as in Sect. 12.11 Since the cluster is considered stationary, 
(8) is no longer relevant. 

The geometry in this case is that shown in Fig. [31 



3. Calculation and results 

3.1. One-dimensional analysis 

For simplicity, we first make a one-dimensional analysis, 
i.e. we consider multiple images only in one direction, effec- 
tively assuming a x spatial hypersurface, hereafter, 
written T"'^. 



3.1.1. Self-gravity of a large cluster 

Figure [H shows the geometry of the situation based on the 
assumptions and general calculation choices of Sect. [2l for 
a 1-torus, model, where, for simplicity, we assume that 
the cluster is moving directly towards one of its adjacent 
images, rather than at an arbitrary direction. A more accu- 
rate calculation would only modify the present calculation 
by less than an order of magnitude. 

The Newtonian attractive force towards the slightly 
closer of the two topological images of the adjacent clus- 
ter, i.e. towards the right in Fig. [2l is then: 



F 



Gm"^ 



1 



(L - x)2 (L + xy 



(2) 



where m is the mass of the cluster, G is the Newtonian 
gravitational constant, L is the comoving size of the fun- 
damental domain in the chosen direction, and a; << L is 



Roukema et al.: Weak Jlj^ in multiply connected universe 



5 



the displacement in the direction of the closest topologi- 
cal image, in physical coordinate units. By the time lapse 
assumption (7), the comoving distance L is 
Let us define 



e := 



F 



X 

L' 
Then 



(3) 



m2 



(1-6)2 (1+^)2 



G 



[4e- 



2e + 3e^ 



.)-(l-2e + 3e2 + ...)] 



L2 



AG- 



(4) 



If we rewrite this as an acceleration and substitute back 
the definition of e, then we have 



4Gr 



L3 



The solution to this equation is the exponential: 



(5) 



(6) 



This is qualitatively what is expected from a cosmolog- 
ical constant: an exponentially growing length scale. 

Could this have any relation to exponential growth in 
the scale factor a(i), i.e. could it provide a cosmological 
constant, or at least a form of dark energy consistent with 
that observed? 

Before discussing this question in the next section, we 
first note that, given the causal contact assumption (4), 
we could expect that not only the closest topological image 
would have an effect on a cluster, but also successive images. 

Eq. ^ for the first N successive pairs of topological 
images then becomes: 



F 



Gm^ 
+ 



1 



1 



(L-x)^ {L + xY 
1 1 



4.8G 



(2L - xf (2L + xf 
2 N ^ 



+ 



AGm 



1=1 



9 

me 
~1X 



(7) 



for N » 1. 

This is only a small correction to Eq. (U) — because 
an effect weakening with the cube of the distance decreases 
rapidly. 

3.1.2. Effect relative to the cosmic web 

As mentioned above, in Sect. 12.21 the geometry for an ob- 
ject slowly starting to fall towards a massive cluster, i.e. 
falling towards a dense node of the cosmic web of density 
perturbations, is that shown in Fig. [3l 



The equation for the long-distance component of accel- 
eration is algebraically the same as for cluster self-gravity, 
i.e. as in Eq. ([2]), except that there is also an acceleration 
term —G^ caused by the local copy of the cluster near 
the test object, since we are interested in long-distance ef- 
fects, and for convenience, we divide by the mass of the test 
object: 



-G- 



-G 



x^ 
m 



Gm 

AGm 
~L3~ 



1 



1 



(L-x)^ {L + xf 



(8) 



The first term in the second line of Eq. ([S]) represents 
local attraction, i.e. Newtonian gravity as it is normally 
thought of, inversely proportional to distance, but the sec- 
ond term is a long-distance term, identical to Eq. ([5]), di- 
rectly proportional to distance. 

However, although this long-distance term is 0/176- 
braically identical to the right-hand side of Eq. ([5]), the 
interpretation is different. 

Instead of the equation representing a high mass cluster 
which has been perturbed from its position and is exponen- 
tially accelerated away from its initial position, in this case 
we have a test particle which feels (in addition to local 
acceleration towards the cluster potential well) an acceler- 
ation away from its nearby (but multiply imaged) cluster 
potential well. 

Again, this is qualitatively similar to the effect of a cos- 
mological constant, since it is repulsive. Moreover, it is an 
effect additional to local gravitational terms — which is 
what would be required of something providing a cosmo- 
logical constant. 

However, since we make assumptions about the cosmo- 
logical epoch, the effect is unlikely to be constant with time, 
so it is more accurate to say that the effect is qualitatively 
similar to the effect of a positive dark energy term. 

Note that, in this case, if we consider only the period 
of infall, before any path crossing or virialisation occurs, 
i.e. when x decreases with time, then force F slows down 
the rate at which the test particle falls towards the clus- 
ter. Nevertheless, the particle is starting to fall towards 
the cluster — in comoving coordinates — so, if an indi- 
vidual test particle has started falling towards the cluster 
in physical coordinates, then the effect from the topological 
images will exponentially decrease with time, as the parti- 
cle approaches the cluster. This is not a problem, since we 
are not interested in following the path of any individual 
particle over time, and since we are most interested in a 
phenomenon which can be related to the system of comov- 
ing coordinates itself, i.e. to the acceleration equation of 
the FLRW model. 



3.2. Three-dimensional analysis 

We generalise the above calculation to an arbitrary dis- 
placement of the test particle from the cluster, in a direction 
not necessarily aligned with the generator g^, and includ- 
ing Newtonian approximation gravitational attraction from 
multiple images generated by orthogonal generators gy and 
gz, for the case of the three-torus T"^. We generalise the 
length L of gx by writing lengths La = L, L^, (Lehoucq 
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Fig. 4. Effect relative to the cosmic web: three-dimensional 
version of Fig. O The small test object is separated from the 
cluster by x along the direction of the generator g^c, and by y 
and z in orthogonal directions. The z direction is not shown. 
Topological images of the cluster ati = ±l,j = fc = are 
shown. 



= L-2[(l + 2e + ...) + (-l + 2e +...)] 

= 4eL-2 

= 4xL-3. (11) 

Unsurprisingly, to first order, this is the same result as in 
Eq. ([5]) . The y and z displacements do not affect the residual 
force in the x direction, they only have second order effects. 

3.2.2. y, z components for i = ±1, j = fc = 

The y and z components of the i ~ ±l,i = fc = terms 
are calculated as follows. 



et al. 1996) for the lengths of gx, gy and respectively. 
Define the displacement vector 



{x,y,z) 



(9) 



of modulus r. 

We continue to assume x,y, z <C La, L^, L„ in order to 
make first order approximations from Taylor expansions. 
The generalisation of Eq. ([5|) is then 

= + ""topo 

= -Gm^ + 



Gm 



E 



{iLa - X, jLe - y, kLu - z) 



(»j,fc)#(o,o,o) - xf + (jXe - yf + {kLu ~ z)'- 



13/2' 



(10) 



where we write i^topo ^^r the residual acceleration due to 
the topological images of the cluster. 

Since the forces add vectorially by orthogonal compo- 
nents, we can first calculate the sum for the same two im- 
ages we used before, i.e. for the two images of the cluster 
along the a;-axis, at i = ±l,j = fc = 0. This is illustrated 
in Fig. m 

3.2.1. X component for i = ±l,j = k = 

The X component of the two terms i = ±l,j = fc = from 
Eq. [iniis, using e := x/La as above. 



^ ((rtopo).=±i. 



j=k=0 



La ~ X 



La X 



[{La - Xf + y2 + z2]3/2 _ ^ ^2 + ^2]3/2 

1 3/2 



3/2 



= L 



(l-e)(l-2e + ...)-^/' 
+ (-l-e)(l + 2e+...)-3/2 

[(1 - e)(l + 3e + ...) + (-1 - e)(l - 3e + 



^ ((rtopo).=±i,.=fe=o)^ 



[(La - XY 2/2 + ^2]3/2 [(_^^ _ ^.)2 + y2 + ^2] 



3/2 



fl-e) 



La 



z 

Ta 



-3/2 



(-1-.)^ + 



2 / \ 2 

y \ I z 



La 



La 



-3/2' 



^yLa 



(1 - 2e + ...y'^'^ + (1 + 2e + ...y^'^ 



-yL-^ [{I + ^e+..:) + {l-2,e + ...)] 



and similarly, 
1 



^ ((i'topo)'=±ij=fc=c 



-2zL~^. 



(12) 



(13) 



Unsurprisingly, this yields a weak net force pulling the 
test object back towards the x-axis joining the two topo- 
logical images of the cluster. 

However, what may be surprising is that this force in- 
creases in amplitude as the test object's y (or z) separation 
increases. For the x-component, it is obvious that as the 
test object becomes less and less symmetrically placed be- 
tween the two topological images, i.e. as x increases, the 
residual force towards the closer image should increase in 
amplitude. But how is it possible that as the test object 
moves further from the x-axis, i.e. as y increases, the resid- 
ual force pulling it back to the plane increases? 

The explanation is in the vectorial nature of the ad- 
dition of forces. (A Euclidean covering space is presumed 
throughout.) The force in the —y direction is only a com- 
ponent of a total, vectorial force. For very small y, the two 
force vectors towards the two topological images almost 
completely cancel since they are nearly perfectly parallel. 
As y increases, these two vectors become less parallel and 
cancel less completely, so although their individual (scalar) 
amplitudes decrease, the y component of their vector sum 
increases. 



3.2.3. Vectorial residual force for x axis and for all images 

Eqs. (|lip . (|12p and (|13p yield the total vectorial residual 
acceleration for the i = ±l,j = k = images of the cluster. 
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i.e. for the closest x-axis images: 



(j^topo)*=±ij=fe=o = Gr 



(4x, -2y, -2z) 



(14) 



By symmetry, the resultant residual acceleration for the 
closest images from all three axes is: 



('^topo)(i,i,fc)6{(±i,o,o),(o,±i,o),(o,o,±i)} 
Gm 



(4x,-2y,-2z]^ (-2x,4y,-2z) 



LI 



2Gm 



LI LI LI 



LI 

(-2x,-2y,4z) 
LI 

2 1 1 

2 11" 
Lf.'a^Ll 



V 



(15) 



If La = Le = Lu, these terms cancel and the resultant 
residual acceleration is zero. 

By symmetry, each successively distant orthogonal, 
equidistant 8-tuplet of topological images in Eq. (fTU]) also 
contributes a zero sum if La = L^ = L„. Numerical calcula- 
tion of the contributions of other symmetrical n-tuplcts of 
topological images distant from the (0, 0, 0) image indicates 
that these also contribute zero to the sum, so that the full 
(first order) sum is zero. 

This shows a three-dimensional effect different to that 
from the x-axis calculation alone: the residual gravitational 
acceleration induced by multiple images disappears if space 
is an "isotropic T'^" model, in the sense that the three 
lengths of the fundamental domain are equal. 

On other hand, in a "slightly anisotropic T'^" model, 
the residual gravity due to multiple images docs not to- 
tally disappear. In order to consider the case in which the 
three side lengths of the fundamental domain are slightly 
unequal, define 



(16) 



Then we have 
2 1 1 _ 

lI'lI'lI 



and similarly 
2 1 1 

Ll Ll L^ 



L-^ [2 - (1 + Se)-^ - [1 + S^) 



-31 



La' [2 -(l-3Je + ...)-(l-3<5n + ...)] 



3L-%5, + 6u) 



3LaH~25e + Su), 



(17) 
(18) 



so that 



(j^topo)(i,j,fe)e{(±i,o,o),(o,±i,o),(o,o,±i)} 
: 6 (^^ + 6^),y {-2Se + S^) , z [S, - 2,5„)] 

For the case La = L^ (i.e. 5e = 0), this becomes 

('^topo)(»j,fe)e{(±i,o,o),(o,±i,o),(o,o,±i)} 
6 GmSuLa^ {x, y, — 2z) 



(19) 



while for the case Le = Lu (i.e. Se = <5«), this becomes 
('^topo)(ij,fc)e{(±i,o,o),(o,±i,o),(o,o,±i)} 



= 6 GmSuLa ^ {2x, -y, -z) 



(21) 



This differs from the acceleration in Eq. by a factor 
of 3i5«. In other words, if the long fundamental dimensions 
are equal to one another and one-third greater than the 
short dimension, then the linearised (in 6u) estimate of the 
effect summed from the three directions, in Eq. (|2ip . is 
as large as if the two larger dimensions were infinite, as 
is effectively represented by Eq. p^ . Clearly, this implies 
that the approximation is valid only for 5u 1/3. 

In each case, this qualitative behaviour is similar to a 
positive dark energy term in the direction of the shorter fun- 
damental length(s), and a negative dark energy term in the 
direction of the longer fundamental length(s), indicating 
that the effect would tend to equalise the three fundamental 
lengths of a model. 

Numerical evaluation of further terms in the sum on 
the right-hand side of Eq. (fTO|) indicates that these modify 
the total sum only slightly. 

4. Discussion 

Is there any relation between the effect found here and dark 
energy? Would this effect really tend to equalise the three 
fundamental lengths of a T'^ model? 

For simplicity, consider the x-axis case and the cluster 
or test object displaced along the x-axis, for the model. 

What has been shown so far is that, under the assump- 
tions listed above, the gravitational effect due to multiple 
topological imaging provides an acceleration proportional 
to displacement, i.e. constant x/x, where either (Sect. [TT]) 
this is the total acceleration for a given massive object to- 
wards its closer topological image, or, (Sect. 12. 2[) it is the 
total long-distance induced acceleration (in addition to lo- 
cally induced acceleration) for a test object "falling" to- 
wards a given massive object. 

The calculation itself is made in comoving space: the re- 
sult of Sect. 12. 1[ for self-gravity of a cluster towards itself, is 
that a perturbed, large massive object is slightly accelerated 
in the direction of its closest image, and this acceleration 
is proportional to the displacement from the initial posi- 
tion in comoving space, so that the displacement increases 
exponentially. Since we are working within the comoving 
frame, it is not obvious, in the Newtonian approximation, 
how to relate this to a modification in the equations for the 
growth of the scale factor itself with cosmological time. 

On the other hand, since the cosmic web is, on average, 
fixed in the comoving frame, it may be possible to interpret 
the second case (Sect. [121), of test objects "faUing" towards 
the most massive objects in the cosmic web, in terms of a 
dark energy term. 

For objects still distant from and falling into massive 
clusters, could the additional force term of Eq. (O, as shown 
in Fig. [3l provide the pressure term in the FLKW acceler- 
ation equation. 



a 
a 



— -p 1 + 3 — 
3 \ pc^ 



(22) 



(20) where a is the scale factor, in order to mimic dark energy? 
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4.1. Amplitude of cosmo-topologica I gravity relative to the 
cosmic web 

An heuristic Newtonian derivation of the acceleration equa- 
tion in this case, for simphcity, treating the x-axis case, for 
the model, follows from Eq. ([8]), where we define x to 
be a fixed length in the comoving reference frame so that 
X = ax- L has already been defined to be a comoving length, 
but in the above equations we implicitly used a = 1; here, 
since we want expressions valid at arbitrary values of the 
scale factor, not only a = 1, we write aL rather than L. 



ax 



X 



X'' 



G 



Am 



Anp Anp Aa^x^ 

'G— ax + G—-^ax 



(23) 



if we estimate that the cluster mass was obtained from mat- 
ter spread at the mean density p throughout a sphere of 
radius x. 

Note that assumptions (5) and (6) are crucial here, since 
we assume that the relevant matter density, both locally 
and at long distance, is that contained inside of local and 
topologically imaged, distant spheres, fixed within comov- 
ing coordinates, around the cluster centre. 

Dividing both sides of Eq. ([23]) by a: = ax yields 



a 
a 



4^ 
L3 



(24) 
(25) 



Using the standard notation for a dark energy compo- 
nent, w :— p/(pc^), we can rewrite this 



4 /X 
3 \L 



(26) 



Since we are interested in objects "falling" (in comoving 
coordinates) towards the nodes of the cosmic web, i.e. at 
most a few tens of Mpc from those nodes, then for a universe 
side length as large as the diameter of the surface of last 
scattering, L — 20h~^ Gpc, we have 



(27) 



So, for the x-axis case in a model, while cosmo- 
topological gravity has the right algebraic characteristics in 
this heuristic approach, with the assumptions listed above, 
to provide an acceleration similar to a dark energy term, 
its amplitude in the present-day Universe is certainly too 
small to be significant, except possibly for extremely high- 
resolution iV-body simulations of the formation of structure 
in the Universe. 

Of course, if the length scale of large scale structure were 
nearly as large as that of the Universe itself, i.e. \i x ^ L, 
then the amplitude of this effect would be much larger. 
However, this would not be physically realistic according 
to our understanding of structure formation. 

As shown in Sect. 13.21 if we consider a displacement of 
the test particle in an arbitrary direction relative to the 
generator vector, and if we sum the contributions from the 
three directions in a T"^ model of nearly equal fundamen- 
tal lengths in the three directions, then a similar heuristic 
argument, using an anistropic scale factor 



implies a similar, but anisotropic, effective acceleration 
term approximately 6 times weaker than for a model 
[Eqs (Uni), ((2T|) ]. where 5 is the fractional difference in 
fundamental lengths. 

In contrast to the T"'^ case, if the three lengths are ex- 
actly equal, then this effect cancels to zero. The fundamen- 
tal lengths must be slightly unequal in order for there to 
be an effect. 



Interestingly, as noted in Sect. 13.2.31 this effect will be 
anisotropic in such a way as to oppose the anisotropy of the 
three fundamental lengths, tending to push the three funda- 
mental lengths towards equality. This is probably the first 
time that such an effect tending to induce La = = for 
a T"^ model has been found. The latter is often assumed for 
simplicity and aesthetic reasons, but here we seem to have a 
physical motivation for this equality as a stable equilibrium 
rather than as an arbitrary assumption. 



4.2. Caveat: Time-varying mass of cluster 

A minor caveat to note for cosmo-topological gravity is that 
the mass m is not constant with time in our model: grav- 
itational collapse will continue and successively larger and 
larger objects will form. However, it has already been noted 
above that this effect is linked to the cosmological epoch, 
and is thus likely to vary with time, so this is not necessar- 
ily a strong argument against some role for this effect as a 
dark energy term: the low amplitude of the effect is a much 
greater problem. 




10 

X (Mpc//i) 



100 



Fig. 5. Ratio of the topological acceleration effect (the effect 
relative to the cosmic web) over the normal gravitational infall, 
shown as a function of the distance of a test point to the cluster 
exerting the greatest gravitational pull on it, using a Hubble 
Volume simulation by the Virgo Supercomputing Consortium 
(see Sect. 14. 4|) . This plot shows the case showing the ratios 
of the modulus of the accelerations. The line shows the analytical 
estimate 4 (f )'' [see Eqs and (f26l) ]. 



a(t) = (aa,ae,a„)(t) 



(28) 
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Fig. 6. As per Fig. (5] showing the (negative of the) 

topological-to-normal ratio of the x axis accelerations in the 
case and the function 4 (^)^. All test points were found (as 
expected) to have opposite signs in the x direction topological 
and normal accelerations; hence, the negative was used. 



10"* - 
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X (Mpc//i) 

Fig. 7. As per Fig. [S] showing the ratios of the modulus of 
the accelerations for the T'' case. The lower scattering of points 
show the case with three equal fundamental lengths and the 
upper scattering of points shows the case with Se = 5u = 0.1. 
The line shows the analytical estimate 12Se (-J)^ [see Eqs (|2H) 
and (f26)l ]. 



4.3. Caveat: x not constant within large scale structure unit 

Another minor caveat, in addition to the anisotropy of the 
effect, is that this effect will vary with distance x from the 
nearest big cluster (node of the cosmic web). 

This implies that the spatially averaged value w ^ 
~ (l) ■^ould be needed before comparing a theoretical w 
value with an observed value, such as the presently esti- 
mated value oi w ~ —1. 

Since most cosmological observations relevant to esti- 
mating the parameters of the metric are in practice aver- 
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Fig. 8. As per Fig. [T] showing the (negative of the) ratios 
of the X axis accelerations in the case with equal fundamen- 
tal lengths (lower "x" points) and unequal lengths (upper "x" 
points). Positives of the ratios, i.e. when a test point is attracted 
by the topological residual force in the same direction as the 
local gravitational force, are shown as circles. This latter sit- 
uation only occurs for the case of equal fundamental lengths, 
where the topological force cancels perfectly according to the 
analytical calculation. The line shows 125^ (^)'^- 



aged out over scales larger than that of large scale structure, 
i.e. on scales ^ 100 Mpc, this is not a problem. 

However, it is interesting to see how the ratio of the 
topological acceleration term to the local acceleration term 
towards the dominant nearby cluster varies as a function 
of the length scale towards nearby clusters. This has been 
presented below in Sect. 14.4] 

4.4. Numerical test 

As a numerical check on analytical calculations of this ef- 
fect, the zero redshift output data file for a "Lambda CDM" 
model of the "Hubble Volume" one billion particle simula- 
tion of the present-day distribution of galaxy clusters made 
by the Virgo Supercomputing Consortium (Evrard et al. 
2002) is used here. 

In other words, a present-day synthetic distribution 
of clusters for a model universe with rim ^ 0.3, f^A = 
0.7,^0 = 70 km/s/Mpc and a CDM (cold dark matter) 
initial power spectrum of density perturbations consistent 
with these parameters, is usedQ The data file contains mass 
and positions of "clusters" detected as virialised groups of 
particles by using a friends-of-friends (FOF) algorithm with 
a linking length of 6 = 0.164. Individual particles have 
masses of 2.24 x lO^^/i-^Mo. 

The side length of the cube is 3/i~^Gpc. This length 
scale is about six times smaller than the likely minimum 
value of L, but is about one and a half orders of magnitude 
larger than the length scale of "large scale structure" , i.e. 

^ See: |http7//w ww . mpa- gar ching . mpg . de/Virgo /bubble . html 

the file is downloadable at tEe time oT submis- 
sion as Ihttp: //In-s .net/E3sl or as http://www.mpa- 
garching.mpg.de / Virgo_Data/hubble_cluster / snapshot /Icdm/- 
cluster_fof_lcdm_z0_b0.164.tar.gz. 
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the scales on whieh clusters form voids and walls, so it 
should be good enough to qualitatively test the result of 
analytical calculations. 

For each test point, the vector accelerations as described 
in Eq. ([TU]) are calculated. That is, the standard gravita- 
tional acceleration —Gm^ towards the cluster whose grav- 
itational pull is strongest at the test point and the "topo- 
logical" gravitational component rtopo f'^'' the topological 
images of this cluster are calculated. This is done: 

(i) for the case (2 adjacent topological images), 

(ii) for the T"^ case (8 adjacent topological images) with 
three exactly equal fundamental lengths, and 

(iii) for the T'^ case with y and z axes' fundamental lengths 
a small fraction Se = Su = 0.1 larger than the x axis 
fundamental length. 

From Eqs p4|) and ((2T|) and the derivation leading to 
Eq. (j26|l . we would expect these three ratios (for the x axis 
direction which dominates, without loss of generality) to be 
approximately : 

(i) Hi)', ^ 

(ii) < 4 (f )^, and 

(iii) 126, (if 
respectively. 

Figures [5][8] show that in a realistic simulation on a 
scale just one order of magnitude smaller than the diame- 
ter of the last scattering surface (Evrard et al. 2002), and 
on scales far enough from individual clusters to be partici- 
pating in the linear regime of density perturbation theory 
(i.e. still expanding with the Hubble flow), the ratio of the 
"topological" acceleration to the normal local acceleration 
is consistent with these three expressions. 

Of course, since the size of the simulation is 3/i^^ Gpc 
rather than a more realistic scale of 20h~^ Gpc, the accel- 
eration ratios in these figures should be reduced by approx- 
imately (20/3)'^, i.e. by a factor of about 300, in agreement 
with the estimate presented below in Eq. ((27|) . 

4.5. Assumptions (5), (6), the topological one-body and 
two-body problems 

Another caveat, probably more important, is that the cal- 
culation we have made, based on assumptions (5) and (6) 
(Sect. 12. ip . is equivalent to what we might call the "topo- 
logical one-body and two-body problems" . 

Naively, it may seem that both the Newtonian and 
relativistic versions of gravity exclude gravitational self- 
interaction of a single, massive, point-sized body: a grav- 
itational interaction between N bodies normally requires 
> 2 bodies. 

However, this intuition relies on the implicit assumption 
that space is an infinite, simply connected Euclidean space. 

In a multiply connected space, a single body can be 
thought of — in the covering space — as a set of multiple 
bodies. Although A^ = 1 bodys exists in a true physical 
sense and in the fundamental domain, A^ ^ 1 bodies exist 
in the covering space, which is the simplest space in which 
to calculate gravitational interactions. 

This is why it is possible for a single body to have 
gravitational self-interactions according to the standard 
Newtonian approximation of gravity. Hence, the "one-body 



problem" is a serious dynamical problem, despite its appar- 
ent absurdity. 

The method of calculating these interactions depends 
on how realistic our universe model is. 

The assumptions (5) and (6) permit a relatively sim- 
ple calculation, and are equivalent to assuming a flat (even 
though mostly empty) space which contains just one (or 
two) objects, i.e. our assumptions are equivalent to a very 
simple, though physically unrealistic, model. 

Also, the distant, supra-SLS images of a cluster are not- 
yet collapsed perturbations, so approximating them as sin- 
gle, point-like objects, while ignoring matter outside the 
radius from which they form is somewhat arbitrary. 

Is it possible to make a more precise calculation which 
avoids assumptions (5) and (6)? Numerical calculations 
would be possible in principle, and would provide a good 
foUowup to our present result. 

However, these are unlikely to be easy, and would be 
different from standard A-body simulations because of as- 
sumptions (7) and (8). Standard A^-body simulations use 
various numerical approximation techniques in order to 
make the calculation time short enough to be practical, 
and, de facto make at least one of the following two ass- 
sumptions: 

(i) instead of assumption (4), assume that something like 
inflation has not occurred 

(ii) instead of assumptions (5) and (6), assume that all 
contributions to the local gravitational potential due to 
long-range, supra-SLS distance scale density perturba- 
tions perfectly cancel out. 

Since inflation scenarios are good candidates for provid- 
ing some ingredient of the correct model of the Universe, 
there is a good chance that assumption (i) is incorrect. 

Moreover, it is not clear to what extent assumption (ii) 
is a good approximation. 

Newtonian gravititational attraction between two point 
objects decreases according to the inverse square of the dis- 
tance, i.e. as oc r~^, but the mass in successive spherical 
shells of equal thickness increases as r^: an anisotropy of 
fixed (small) solid angular size cos Odcj) d9 at all radii would 
provide equal attractive forces from each successive shell, 
no matter how distant. 

This is no problem in a perfectly homogeneous universe: 
hence, the FLRW solution of the Einstein field equations. 

But in a universe with perturbations, i.e. a "slightly" 
inhomogeneous, anisotropic universe, the contributions do 
not perfectly cancel. 

In a nearly FLRW universe, i.e. one containing a spec- 
trum of density perturbations, the degree to which the long- 
distance terms cancel most likely depends on the full nature 
of the perturbation spectrum. 

Farrar & Melott (1990) suggest from their Af-body cal- 
culations that assumption (ii) is correct in the case of a mul- 
tiply connected fiat space (of equal fundamental lengths). 
Further work in this direction would be interesting, in par- 
ticular by introducing the physically standard assumption 
(not used explicitly in normal A^-body simulations), i.e. our 
assumption (7), that cgw is finite and of the standard value 
CQW — csT ~ CQR = cem- Testing the effects of slightly 
unequal fundamental lengths would also be interesting. 
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4.5.1. Inhomogeneities and the Einstein-Hilbert equations 

As has been pointed out for many years, e.g. Section 7, Ellis 
& van Elst (1999) and references therein, finding an exact 
inhomogcncous solution to the Einstein-Hilbert equations, 
and then averaging this to find a "mean density" homoge- 
neous solution, is not mathematically equivalent to adding 
perturbations to the exact FLRW (homogeneous) solution 
to the Einstein-Hilbert equations. The latter, approximate 
approach gives a good match to observations, but this docs 
not prove that the difference between the exact and pcrtur- 
bative approaches is negligible. 

Buchert & Carfora (2002, 2003) have been studying this 
problem in depth for several years and have noted that the 
difference between the perturbative approach and the exact 
(but still developing) approach leads to a "backreaction" 
term (and also a curvature term) which may have an ef- 
fect on scales ranging from that of a massive galaxy cluster 
possibly up to super-Hubble length scales. Recent discus- 
sions (Kolb et al. 2005b, 2005a) indicate that this presently 
seems, pending full calculation, to be a viable candidate 
for a dark energy term, without the need for, e.g. adding 
any scalar fields, based on the backreaction term from sub- 
Hubblc length perturbations. 

In a qualitative sense, this is similar to the effect dis- 
cussed in this paper, except that we are primarily interested 
here in the effect of perturbations on gravity in the case of 
non-trivial global topology. 

Buchert & Ehlcrs (1997) and Ehlcrs & Buchert (1997) 
presented work closely related to what is presented here, 
considering both T'^ and K"^ models, for the Newtonian case 
of perfectly isotropic expansion: 

aait) = aeit) = auit) Vi. (29) 

For the case, they found a consistent result with our 
first order result, i.e. that for a compact space admitting 
a global Hubble flow, i.e. if the average shear and rotation 
are zero (Appendix A, Buchert & Ehlers 1997; sec also four 
equivalent statements related to this in Sect. 3.4 of Buchert 
& Ehlers 1997), then the perturbations have no effect on 
the rate of expansion according in the Newtonian approxi- 
mation. 

Given that our heuristic calculation suggests that for a 
slight anisotropy in the three lengths scales of T'^ , the ex- 
pansion rate should be anisotropic, it would be interesting 
to develop Buchert & Ehlers (1997)'s approach further us- 
ing (slightly) anisotropic scale factors, i.e. using Eq. (^5)1 to 
replace Buchert & Ehlers (1997) 's equation V{t) ~: aj){t) 
(just after eq. (3)) by 

V{t) = aaa,au{t) (30) 

and following through, using the definition in Eq. (jl6p . The 
related work in Ehlcrs & Buchert (1997) could be reworked 
starting from their eq. (15). 

Just as perfect homogeneity is certainly wrong, perfect 
isotropy in the expansion rate is probably also just an ap- 
proximation to a more accurate model, whether or not the 
correct model is closer to T'^ or or another model. 

5. Conclusion 

A residual gravitational effect, which wc could possibly call 
"cosmo-topological gravity", occurs due to distant multi- 
ple topological images in a multiply connected universe 



which "remembers" the gravitational potential generated 
by multiple topological images, which in the covering (ap- 
parent) space arc located outside of the present surface of 
last scattering, and were causally contacted at some earlier 
epoch, for example due to some prior amount of inflation. 
A Newtonian approximation, in which the speed of trans- 
mission of gravitational information is flnitc (equal to the 
special relativistic space-time constant c), rather than infi- 
nite as in a fully Newtonian calculation, is used here. 

For a low mass test object "falling" (in comoving coordi- 
nates) towards a relatively nearby (at a comoving distance 
X less than a few tens of Mpc) large, massive collapsed ob- 
ject at the present epoch, i.e. a massive cluster of galaxies, 
in a 3-torus universe of side-length i « 20 Gpc, the 
two closest topological images of the cluster together yield 
a residual Newtonian force on the test particle which locally 
appears as a force repelling the test particle away from (its 
nearby image of) the cluster, provided that we consider 
the test object lying along the axis joining the two closest 
topological images of the cluster in one of the fundamental 
directions, and wc only consider these two images. 

This residual force provides an acceleration alge- 
braically similar to that of dark energy, but weaker by 
many orders of magnitude, i.e. by a factor of approximately 
^ at the present epoch. 

A more general, three-dimensional calculation, for a test 
object displaced in an arbitrary direction in a 3-torus uni- 
verse, shows that the effect cancels out to zero if the three 
side lengths are exactly equal. If the side lengths are slightly 
unequal, by a fraction ~ 6, then an anisotropic dark energy 

term, about S (^) ^ 10~^ times weaker than the observed 
dark energy (for side length L « 20 Gpc), will accel- 
erate the expansion of the shorter length(s) and decelerate 
the expansion of the longer length(s), tending to equalise 
them. 

This is probably the first known physical effect which 
could relate the three fundamental lengths of a universe 
with T"^ spatial sections. The equality of the three lengths 
in simulations of the T"^ model has often been assumed, but 
on purely aesthetic grounds, without any physical justifica- 
tion. 

It is clear that the effect is not significant in the present- 
day Universe. 

Apart from developing an anisotropic version of Buchert 
& Ehlers (1997) 's approach, an interesting (and challeng- 
ing) foUowup project would be to check whether or not 
the same effect occurs in spherical, multiply connected uni- 
verses: could the residual gravitational force due to multiple 
imaging have helped push the shape of the Universe into 
that of the Poincare dodccahcdral space (PDS) during or 
not long after the quantum epoch? Did it help isotropise 
the Universe? 
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